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ABSTRACT 

The Mid-IR colors (Fg/F^) of galaxies together with their IR-UV luminosity 
correlations can be used to get some insight into the relative abundance of the 
different dust grain populations present in them. The ELAIS-N1 field contains 
thousands of galaxies which do not have optical spectra but have been observed 
in the Mid-IR by Spitzer and UV by GALEX making it ideal for these studies. As 
part of this work we have selected a sample of galaxies from the ELAIS-N1 field 
which have photometric observations in the MIR and UV as well as photometric 
redshifts from the SDSS database. We put the constraint that the redshifts are 
< 0.1, thereby giving us a total of 309 galaxies. We find that the majority of the 
galaxies in the sample are PAH dominated due to their high MIR flux ratio. We 
also find a reasonable correlation between the Mid-IR and the UV luminosities out 
of which the Mid-IR emission from PAHs at 8 /mi is marginally better correlated 
than the 24 /mi VSG emission with the UV luminosities. However, if we divide 
the sample based on their F 8 /F 2 4 ratios which is also an indicator of metallicity, 
the MIR-UV correlation seems to increase with the F 8 /F 2 4 ratio. But the MIR- 
UV correlations are not very different for the PAHs and the VSG population 
within the individual metallicity groups. 
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Introduction 


Most of the mid-infrared (MIR) emission (3-40 /irn) from star-forming galaxies is due 
to thermal emission from dust grains heated by stellar radiation. The dust emission consists 
of two different components, i.e., broad emission features at 3.3, 6.2, 7.7, 8.6 and 11.3 
/irn from large aromatic molecules known as polycyclic aromatic hy drocarbons or PAHs 
which are known to have vibrational transitions at these wavelengths (ILeger and Puget 1984 
Alla/mandola et al.l Il985 1 and con tinuum emissi on from very small dust grains (15-40~A, 
sometimes referred to as VSGs) flSellgrenl 11984 ). The PAHs are excited by single photon 
heating and their Mid-IR luminosities are therefore independent of the radiation field density 
(Draind 20031) . The 8 /irn emission from PAHs is correlated with cold dust emission at 16 0 pm 


(Bendo et ah] 20080, making them good tracers of the diffuse ISM in galaxies (Peeters et al. 


2004 b Laboratory studies have shown that PAHs are excited mainly by UV photons (see 
Sellgrenl (120011 ) and references therein) and there have been several observations of these 
emission features in regions of intense UV radiation. However, the detection of these emission 
features in benign environments with little UV radiation as well as in some early-type galaxies 
sho we d that they could also be excited by evolved stellar populations (for a recent review see 
Calzetti (2011)). Such regi ons are d omin ated b y large PAHs (N c > 100) where the absorption 
extends into the visible (ILi and Drainell2002al ). The 8 pm emission from PAHs is found to 
be strong outside HII regions where individual PAH molecul es can s urvive an d t h ey are 


irrad iated by the general interstellar radiation field (ISRFi flGordon et al.ll2008t iHelou et al. 
2004h . 


The 24 pm continuum emission from VSGs (which may be PAH clusters or small car¬ 
bonaceous grains), can only be excited by highly luminous UV em itting stars and are found 


to be associated with the inner parts of HII regions. Zhu et ah (12008! ) have found a tight 


correlation between the 24 pm (continuum emission from VSGs) and 70 pm luminosities 
both of which represent the hot dust. 

The 8 pm emission from PAHs has be en found to be tight l y correlated with metallicity 


of th e environment in which it is present flDraine and Li! 120071: ISmith et al.l 120071: IWu et al. 
20071: iHunt et al.l 120101). This holds both globally and locally as in the case of HII regions 


( Gordon et ah 20081: Iwiebe et al' 2011 : Khramtsova et al. 2013 . 2014 ). where a tight correla¬ 
tion has been found between the F 8 /F 2 4 ratios and the metallicities. While higher metallicity 
galaxies or regions lead to higher PAH abundances, lower metallicity objects, for eg., the 
Small Magellanic Cloud (SMC), tend to have much lower abundances with < 0.4% of the 
interstellar carbon residing in PAHs (ILi and Draine 2002bl ). This may be either due to lower 
rates of formation of PAHs in low metallicity objects or due to higher destruction. Low 
metallicity objects tend to have strong radiation field that easily destroys the PAHs except 
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the very large ones or PAH clusters (IMadden et al.ll2006l: iRapacioli et al.l 120061) . Thus, low 
metallicity environments with associated high radiation held are dominated by emission from 
VSGs compared to PAH emission. Several studies have found a reduction in the F 8 /F 2A ra¬ 
tios in high radiation held regions which is lower than the predicted dust model values and 
have attributed it to the destruction of PAHs in intense radiation held regions or due to an 
increase in the VSG population due to the destruction of larger grains or both. However, 
the F 8 /F 2 4 ratio is also found to be age dependent for different values of metallicities of HII 
complexes ( Khramtsova et ah 2014). Therefore the line and continuum emission observed in 
different interstellcr environments are associated with different compositions of dust. 


Using a combination of the MIR emission from VSGs and PAH molecules together 
with the observed UV hux from stars in the galaxy, we intend to classify galaxies based on 
their F 8 /F 24 ratios and understand the nature of dust present in them, in the absence of 
spectroscopic data. Finally we compare our results with those obtained for other well known 
galaxy samples and HII regions. 


SAMPLE SELECTION 


We have used observations of the SWIRE ELAIS-N1 held ( RA=244°, Dec=+54°) from 
the Infrared Array Camera (IRAC: 3.6, 4 .5, 5 . 8 fe 8 .0 pm) and the Multi-Band Imaging 
Photometer (MIPS: 24 pm) on the Spitzer ( Werner et al.l 2004 ) and the far ultraviolet (FUV: 
1350-1750 A) and near ultraviol et (NUV: 1750-2800 A) imagers on the Galaxy Evolution 
Explorer ( GALEX ; M art in et ah ( 2005 )). At low redshifts, the IRAC 5.8 pm and 8.0 pm 
channels detect PAH emission features, with theJhO pm channel detecting it up to a redshift 
of 0.2 due to its larger bandwidth (IHuang et al.l 20 071 ). The other two IRAC channels 3.6 pm 
and 4.5 pm detect mostly stellar emission while the MIPS 24 pm band captures the thermal 
continuum emission of VSGs. We have selected those extended sources (Source extractor ext 
hag >= 2 , Bertin and Arnouts (|l996n ) that have been observed in all the IRAC channels, 
in the MIPS 24 pm band and the GALEX FUV and NUV bands. There are a total of 1169 
extended sources in the held that have significant emission at the wavelengths considered 
here. We consider only those sources which have photometric redshift < 0.1 as observed 
by SDSS. With this restriction we are left with only 309 extend ed so ur ces. We have used 
the Kron huxes from the SWIRE 2005 catalog for the data (ILonsdale et al.l 12003 ) and the 
GALEX huxes from the GALEX band merged catalogue corresponding to the tiles ELAIS- 
Nl-(00,03,12,01,02). The UV huxes have been corrected for Galacti c ext inction using the 
E(B-V) values from the Spitzer cata l og wh ich originate from lSchlegel et al.l (119981) . We have 


used the co-efhcicnt of Helou et al. (120041 ) to correct for the stellar contamination in the 8 
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fim band. 


3. F$/F ‘24 ratios and Dust model comparison 


We compared the sample of galaxies in the ELAIS-N1 field with SDSS observed sources 
and found 309 of them to have photometric redshifts. We then converted the redshifts to 
distances and computed Luminosities for these sources. We have plotted the L 8 /L 2 4 ratio as 
a function of number of galaxies in Fig. [0 We can see that majority of the galaxies have a 
high ratio indicating a higher PAH content in them. In order to study galaxies with different 
relative abundance of PAHs, we have divided the sample into four groups (R1-R4) based 
on their L 8 /L 2 4 ratios. These groups can be considered to represent different metallicities 
since the F 8 /F 2 4 ratios are found to be tightly correlated with galaxy metallicities. Galaxies 
belonging to different groups are marked by different symbols in Fig. [0 


We have compared the observed F 8 /F 2 4 ratios with the theoretical calculations of 


Draine and Li (120071 ). where they have used a mixture of carbonaceous grains (comprised of 


PAHs and graphitic particles) and amorphous silicates exposed to radiation fields ranging 
from 0.1 — 30000 U, to predict rnid-IR fluxes in each of the instrumental bands of Spitzer. U 
is the local interstellar radiation field in our Galaxy (Math is et al. 1983) • For each radiation 
field model, the intensities are also computed for different values of the PAH mass present 
in grains denoted by qPAH, which is a function of the metallicity. 


We have plotted the median of the observed L 8 /L 24: ratio vs. L 8+24: / Ljjv ratio for each 
group (denoted by ”x”) together with the scaled model values, where U is taken in place of 
the UV fluxes in Figures [2] and El The same values of U have been used for the FUV and 
NUV fluxes in order to compare the model with the observations. The model values for the 
two qPAH fractions 4.58% and 0.47% are marked as squares and triangles respectively. Even 
though the local interstellar radiation field density, U has a spectrum which has a peak in 
the visible due to all the late type stars in the Solar neighbourhood, we can approximate 
it by the UV fluxes when studying the MIR flux ratios because the dust grains and PAHs 
are known to have very high extinction cross-sections in the UV compared to the optical. 
We see a reasonable fit between the model and observed profiles especially for the high PAH 
content. This may be an indication that the dust model with a high PAH fraction of 4.58% 
is a fairly good representation of the dust content of these galaxies. Note that this value is 
representative of Milky Way dust as well. 
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Fig. 1.— Four groups of galaxies are represented by different symbols in the figure. The 
classification is done based on their L 8 /L 2 4 ratio. 



/ L FUV 


Fig. 2.— Logarithmic plot of the median of observed log(L 8/i /L 2 4 M ) vs. log(L 8M +L 2 4 ^)/-b,Fi/v 
(cross) and model flux ratios with U taken in place of FUV. The model values for qPAH = 
4.58% and qPAH = 0.47% are plotted as boxes and triangles respectively. 
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4. CORRELATION STUDIES 


As mentioned earlier, for 309 galaxies of our sample we have converted the photomet¬ 
ric redshifts to distances and computed Luminosities. We studied the UV-IR luminosity 
correlations of our sample and looked for variations with metallicity. 

We find a good correlation between the UV and the 8.0 /irn IRAC luminosity (see Fig. 
2](a) and (b)), which is most sensitive to the PAH emission, with a better correlation for the 
NUV luminosities. Much of the scatter seen in the figures is probably due to differences in 
the intrinsic UV extinction, which will depend on the dust content, geometry and metallicity 
for each galaxy and hence difficult to correct for in this sample due to the unavailability of 
spectra. In Fig. [4] ((c) and (d)) we have also shown the correlation between the UV and 
8.0 /irn IRAC flux. The 24 /irn is also correlated with the UV luminosities as seen in Figure 
[5] The Spearman Rank Correlation co-efficients for these sources are tabulated in Table 
m We see that the 8 /irn emission is marginally better correlated with the FUV and NUV 
luminosities than 24 /irn. This together with the higher F 8 /F 2 ^ ratio for the majority of the 
galaxies shows that the PAH emission dominates in these galaxies and the PAHs are mostly 
excited by the UV emitting stars present in them. Since the 24 /irn feature requires VSGs 
to be heated by very hot massive stars, the slightly lower correlation together with the high 
F 8 /F 2 4 ratios could indicate that in this sample, most of the galaxies do not have a large 
number of hot massive stars or have high metallicities as also seen from the comparison with 
the dust model. 


In order to verify if the correlation co-efficients vary with the metallicity, we now compute 
the Spearman Rank correlations for the five 8 fim /24 fxm groups seperately. We see that the 
correlation co-efficients reduce marginally with decrease in MIR flux ratios or metallicities 
i.e. as we go from group 1 to 4. But there is no difference in correlation between the PAH 
and UV emission or VSG and UV radiation. This is in agreement with studies of starforming 
galaxies where the correlation co-efficient s between MIR and UV luminosities are similar for 
PAH and VSG emission f Zhu et ah 2008 b 


Table 1; Spearman-Rank Correlation co-efficients between FUV and MIR luminosities for the SDSS sample together with 


the probabilities. Correlations are also given for the different groups seperately. 


Bandl 

Band2 

Rtot 

Probt ot 

Ri 

Probi 

r 2 

Prob 2 

Rs 

Prob 3 

R 4 

Prob 4 

FUV 

8.0 

0.57 

3.2x 10 -28 

0.60 

1.5xl0 -9 

0.57 

l.OxKT 7 

0.55 

2.7xl0 -8 

0.44 

2 . 0 xl 0“ 4 

FUV 

24.0 

0.55 

5.4xl0 - 26 

0.61 

8 . 6 xl 0 _1 ° 

0.57 

1 . 2 xl 0 -7 

0.55 

2.3x10-® 

0.42 

4.3xl0 “ 4 

NUV 

8.0 

0.64 

9.4x 10 -37 

0.66 

8.5xl0 -12 

0.64 

4.3xl0 _1 ° 

0.62 

l.OxlO - 10 

0.54 

4.2X10 " 6 

NUV 

24.0 

0.62 

1 . 8 x 10 -34 

0.66 

4.4xl0 -12 

0.64 

5.2xl0 1 ° 

0.63 

8.4X10 - 11 

0.51 

l.lxlO -5 
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5. Conclusions 


We have used archival Spitzer and GALEX data of the integrated light of galaxies in the 
SWIRE ELAIS-N1 held in order to study and classify these galaxies on the basis of the type 
of dust grains in them. We have used the F 8 /F 2 4 ratio to deduce the relative abundance 
of the PAH and VSG population. We find that majority of galaxies in the sample are 
PAH dominated. Correlation studies of these galaxies show that that the 8 /jrn emission is 
reasonably correlated with the FUV and NUV radiation. This is also true for the 24 /jrn 
emission. By considering the observed F 8 /F 24 ratio as a proxy for the metallicity of the 
galaxy, we have divided the sample into four groups of different metallicities. We find that 
the observed MIR flux ratio follows the dust model ratio reasonably well especially for the 
higher PAH fraction. 


Depending on the dominating dust component (PAH or VSG) that absorbs/scatters the 
UV radiation we could see a change in the MIR-UV correlations. To check this, We have 
compared the MIR-UV correlations for the high and low metallicity galaxies, and find that 
the correlation co-efficients between the different groups reduce with reduction in the MIR 
flux ratio. For low metallicity galaxies, the harder radiation field flMadden et al.l 20060. could 
be responsible for destruction of the PAH and VSGs thereby affecting the correlations. Or 
this could be just because of the higher extinction in these galaxies due to the larger dust 
content. However, the correlation co-efficients are not very different for the PAH and VSG 
emission in each group, thereby making it difficult to conclude which grain population affects 
the UV radiation field in these galaxies. This could be because the abundance of the different 
grain populations is less sens iti ve to t he UV radiation field than the metallicity. This has 
been seen in HII regions (Khramfsova et al. 2014) and has been attributed to production of 
PAHs by aromatization of carbon grains in regions of high radiation field or low metallicities. 
Similarly the VSG population could also increase due to destruction of larger grains in these 
regions thereby influencing the MIR flux ratios. 


Hence using this method we can see that even in the absence of spectroscopic mea¬ 
surements we can study the relative abundances of the different populations of grains in 
galaxies and the results are consistent with those obtained for well studied galaxies and HII 
regions with metallicity measurements from spectra. However, it is not possible to identify 
what factor alters the relative abundance of the different grain populations in these galaxies, 
which probably means that several factors influence their relative abundance both globally 
and locally. 


This work is based in part on observations made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology under a con- 








tract with NASA. This article is also based on data obtained by the NASA mission Galaxy 
Evolution Explorer (GALEX). PS would like to thank DST Fast Track programme (grant no.: 
SR/FTP/PS-68/2008) for financial support. RG would like to thank the IUCAA Associate- 
ship programme for their support and hospitality. AP thanks DST-SERB FAST TRACK 
programme and DST-JSPS project for financial support and IUCAA for visiting associate- 
ship. 
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Fig. 3.— Logarithmic plot of the median of observed log(L 8At /L 24At ) vs. log(L 8/i +L 2 4 M )/Ljv; 7 v 
flux ratios (cross) and model flux ratios with U taken in place of NUV. The model values 
for qPAH = 4.58% and qPAH = 0.47% are plotted as boxes and triangles respectively. 
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Fig. 4.— (a) and (b) represent the logarithmic plot of the observed UV vs 8 /jrn luminosities 
and (c) and (d) represent the logarithmic plot of the observed UV vs 8 //m fluxes. 
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Fig. 5.— (a) and (b) represent the logarithmic plot of the observed UV vs 24 /mi luminosities 
and (c) and (d) represent the logarithmic plot of the observed UV vs 24 /mi fluxes. 







